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Abstract. - We show that an independent four-body momentum scale /i(4) drives the tetramer 
binding energy for fixed trimer energy (or three-body scale ^(3)) an d large scattering length 
(a). The three- and four-body forces from the one-channel reduction of the atomic interaction 
near a Feshbach resonance disentangle fj,^ and /U(3). The four-body independent scale is also 
manifested through a family of Tjon-lines, with slope given by M(4)/^(3) for a -1 = 0. There is 
the possibility of a new renormalization group limit cycle due to the new scale. 



Recent progress in experimental techniques, creating tunable few-body interactions in 
trapped ultracold gases [1], is allowing the investigation of few-body quantum phenomena, 
with long wave-length triatomic properties being currently probed. In recent experiments [2,3] 
it was confirmed the universal properties coming from the Efimov physics [4] occurring for 
large two-body scattering lengths [5-9] . So far experimentally, it was discovered that magnetic 
tunable molecular interactions in traps allow the formation of complex molecules like Cesium 
tetramers (Cs 4 ) near the scattering threshold [10]. One can foreseen that the physics of 
four-atom systems in ultracold gases will be probed in the future with tunable interactions. 

Near a Feshbach resonance the two-atom scattering length, a, can vary from very large 
negative values to positive values, allowing virtual or weakly-bound dimers. The scattering 
length is large in respect to the atom-atom interaction range (ro), driving to the use of 
concepts developed for short range interactions and halo states [11]. In the limit of large a, 
the interaction can be taken as of zero range [12]. The appearance of Thomas-Efimov states 
in three-boson systems is controlled by the ratio ?"o/|a| ~~ > [4, 13-15]. In this exact limit, 
it is observed an infinite sequence of three-body bound states [16,17], identified [15,18] with 
an underlying renormalization group limit cycle [19]. The collapse of the three-boson system 
when the two-body interaction range goes to zero for fixed a demands one three-body scale 
to stabilize the system. 
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In the nuclear physics context, it was found that the 4 He energy is determined by the triton 
one, which is known as the Tjon-line [20]. Once the triton energy is fitted to its experimental 
value the Tjon-line gives the observed 4 He binding. That result could be particular to the 
nuclear potentials models used in the calculations and may not be valid near a Feshbach 
resonance in atomic systems. 

In this letter we show that a new parameter/scale governs the properties of four- atom 
systems near a Feshbach resonance. We also discuss that the three and four-body energies 
can be disentangled due to the presence of few-body forces from the one-channel reduction. 
First, we discuss the atomic interaction near the Feshbach resonance and its one-channel 
reduction. Second, we solve the four-boson problem for a zero-range interaction with two and 
three-boson binding energies fixed to investigate the necessity of a new scale to determine the 
four-boson binding energy. 

The accepted wisdom is that near a Feshbach resonance only the two-body scale is moved 
through the large variations in the atom-atom scattering length. However, surprisingly enough, 
there are three and four-body potentials induced close to a Feshbach resonance, which can 
disentangle the corresponding scales and justify a thoroughly discussion of the independence 
of these scales for trapped atoms. To understand the origin of the induced interactions, we 
start with the Feshbach decomposition in P and Q spaces (P + Q = 1), with the Q— space 
representing the pair in the potential well where it is bound, while the P-space represents the 
pair in the lower potential well where the continuum channel of the ultracold regime happens. 
Three- and four-body effective potentials appear when the pair of particles forms the bound 
state in the Q— space when interacting with spectators. To make it concrete, we chose an 
example for a three-body system (ijk) where a pair (ij) forms a bound state in the well corre- 
sponding to the Q-space. The Hamiltonian is H = PHP + QHQ + Pv l3 Q + QvjjP, where the 
transition potential between P and Q spaces is v. To simplify the discussion it is assumed that 
the channels Hamiltonians are PHP = H + V ik + Vjk + V^ and QHQ = H + V ik + Vjk + V® . 

H a is the kinetic energy operator. The gap AE = Hindoo (v®(r) — Vij(rfj > gives the dif- 
ference in the asymptotic values of two potential wells. The potential Vij is the non-resonant 
part of the interaction between i and j, while V® is the second potential well where the 
Feshbach resonance lies. The effective Hamiltonian acting in the P-componcnt of the wave 
function is 



where the resolvent Gqq(E) — [E — H — V® — Vik — Vjk] -1 - The last term in the effective 
Hamiltonian carries the resonant term of the potential between particles (ij) and the induced 
three-body potential, as can be understood from the decomposition of the Q— channel resolvent 
as Gqq(E) = [E - H a + V?]- 1 [1 + (V lk + Vjk)G QQ (E)]. The resonant interaction between 

particles ij is given by PvijQ[E — Hq — V®]~ 1 Qv\jP ', when it is accounted the dominance of 

the projection over the bound state of the potential V® near the scattering threshold. The 
three-body connected part of the remaining term corresponds to a three-body interaction, 
with intensity clearly depending on the position of the Feshbach resonance. Such, reasonings 
are easily extended to four particles. Therefore, in principle not only the scattering length 
can be tuned near a Feshbach resonance but also the three and four-body binding energies. 
Turning off the excitation of the Feshbach resonance in the presence of spectator particles, 
the formalism is reduced to the one discussed in Ref. [21]. 

In short, near a Feshbach resonance the induced one-channel three and four boson forces 




(1) 



M. T. YAMASHITA, LAURO TOMIO, A. DELFINO and T. FrEDERICO: FOUR-BOSON SCALE NEAR A FESHBACH RESONANCE3 



can drive independently the corresponding physical scales. Therefore four-body observables, 
like four-boson recombination rates or atom-trimer or dimer-dimer scattering lengths can 
exhibit correlations not constrained by one low-energy s-wave three-boson observable and a, 
if a four-boson scale really exists. 

We begin our discussion of the relevant scales in the few-boson system reminding the form 
of the regulated trimer bound-state integral equation, in units of h = m — 1 (m is the boson 
mass), which is written as 



/(<?) 
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where the second term in Eq. (J2J) brings the physical scale, (j,r 3 ) ~ r^ 1 , to the three-boson 
system [22] producing a finite ground-state energy and avoiding the Thomas collapse [13]. 

The sensitivity of three-boson s-wave observables to the short-range part of the interaction 
in weakly bound systems is parameterized through the value of the trimer binding energy which 
corresponds to the regularization scale /i( 3 ) . The s-wave three-boson observables are strongly 
correlated to the trimer energy and scales in the general universal form [16,23]: 



,(2) 
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where 3 can represent an observable at an energy E or an excited trimer energy (the depen- 
dence on E does not appear in this last case). The exponent -q gives the correct dimension 
to 3 . Eq. J2J is renormalization group (RG) invariant with its kernel being a solution of a 
non-relativistic Callan-Symanzik differential equation [24] as function of a sliding /j,r 3 ) . In that 
way E 3 and three-body observables are independent of the subtraction point (see Ref. [23] for 
a discussion of the RG invariance in three-body systems). 

In a four-boson system with short-range interactions one may be tempted to think in an 
independent momentum scale, which is still under discussion [25,26]. In the hyperspherical 
expansion method used in Ref. [11,25] it is enough to fix the scattering length, effective range 
and shape parameter in the limit of zero-range interaction to stabilizes the N-boson system. 
On the other hand within effective field framework it is enough a repulsive three-body force 
that stabilizes the trimer to allow finite values of the tetramer binding energy [26] . 

In our strategy to verify the necessity of a new four-boson scale, the dimer binding energy 
is set to zero while the trimer ground-state energy is kept fixed. Within these assumptions, 
to allow finite results for the tetramer energy, the Faddeev-Yakubovsky (FY) equations [27] 
should be regularized with the insertion of momentum cutoffs or subtraction scales. From 
our results for the binding energies in the limit of a zero-range interaction, it emerges the 
necessity of a new scale in the four-boson system, which depends crucially on the four-body 
scale, the subtraction point [i 2 ^ . The existence of an independent four body scale can be in 
principle checked by measuring s— wave observables for three and four atoms near a Feshbach 
resonance. If the three and four-body scales can be made truly independent due to multi-atom 
forces, one can distinguish from the case that such scales are equal considering the correlation 
between two observables. 

In the next, we show the systematical inclusion of the three and four-body subtraction 
points in the unregulated FY equations. First, it is natural to require /X( 3 j as the subtraction 
scale of the embedded three-boson system in the four-body one, otherwise the threshold for 
the four-boson energy would be different from the ground-state trimer energy. Second, the 



4 



EUROPHYSICS LETTERS 




Fig. 1 - Jacobi coordinates for the Faddeev-Yakubovsky amplitudes. K-type, left-side, and //-type, 
right side. 

four-boson scale or the new subtraction point is introduced in the set of FY equations, and 
the tetramer binding energy is calculated for different values of /j,u\ . 

The dynamics of a four-boson system can be described in terms of the Faddeev-Yakubovsky 
amplitudes [27], K and H, represented schematically in Fig. Q). The amplitude k corre- 
sponds to the partition in a three-body subsystem formed by particles, ijk, and the spectator 
one, I, while HijM describes the partition in two subsystems formed by the pairs ij and kl. 
There are 18 independent FY amplitudes, 12 if-type and 6 -H-type. The bound state FY 
amplitudes satisfy: 

\K\ hk ) = G Q t ij {E-E ij>k -E l )[\K\ k J + \K] k J 

+ 1*4,*) + \K) k ,) + \H ikJl ) + \H jktil )] , (4) 
\H im ) = Got^E - E im - E kl )[\K{ Ul ) + \Kl u ) + \H kl<ij )] , (5) 

where the four-body free resolvent is Go — (E — i?o) _1 (-^o is the kinetic energy operator) 
and tij is the two-body T-matrix. The kinetic energies for the Jacobi momenta for identical 
particles of unit mass are E ijtk = fg^fc, Ei = E ijtk i = \q% M E kl = ill- Tne total 
wave function of the four-body system is given by the sum \^) — Yl^j where the 

Faddeev components = GoVij\*I>) can be written in terms of the FY amplitudes as 

l^i) = l^i, fc ) + + \H ijM ). 

The renormalized two-body T-matrix for a zero-range interaction is given by tij(x) — 
\Xij) T ( x )(Xij I w ith r_1 ( a; ) = 27T 2 [a -1 — \/—x\ and (qij\xij) = 1- The structure of the FY 
equations is simplified for the zero-range interaction once the separation \K\^ k ) = Go\xij)\fc\j k ) 
and \Hij t ki) = GQ\xij)\Hij,ki) 1S recognized in Eqs.QJ and ©. The FY reduced amplitudes, K 
and 7i, satisfy a coupled set of integral equations which needs regularization. One recognizes 
that the resolvent of the immersed three-boson subsystem should carry the regularization scale 
/i( 3 ) to avoid the Thomas collapse of the system. This is done in analogy with Eq. @, where 
the momentum integration is regularized by a subtracted Green's function in a subtraction 
point /i( 3 ). In the l.h.s. of Eq. © one sees the correspondence with the three-boson equa- 
tion written as T- l (e { l' Z) )\K,\ jk ) - Q^] ik \lC\ kJ ) - g$ jk \K. l jk J = 0. However, other terms are 
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present in the four-body equations which require regularization. 

We introduce the scale /i( 4 ) such that the four-body free Green's function Go (£4) is sub- 
stituted by Ga{E/j) — Go{— M(4)) m a direct generalization of Eq. (J2J), as suggested by [22]: 

G$ ik Wik,i) + IWiM>] + sSl„ [\K* Jktl ) + \n jk , u )} , (6) 

T-H£ ,2) )\Hiw) = G\%t [\K,y + \lC khj ) + \H kl>ij )] , (7) 

(13) (2^) 

where e 2 ' — J5 — J5« ^ — E\ and e 2 = E—E^^i — Em are the two-body subsystem energies 
in the (1,3) and (2,2) partitions, respectively. The projected Green's function operators are 

Qif-ik := (Xij\G^\Xik) w hh N equal 3 or 4, with the subtracted Green's functions given by 

G 3) = [E-H ]- 1 - [-nf 3) -Ha]- 1 and G 4) = [E - H ] - 1 - [-/x^ 4) - Hq}- 1 . Maintaining the 
two- and the three-body scales fixed, the four-body scale can be moved and its consequences 
for the tetramer binding investigated. 

The four-boson integral equations for the reduced FY amplitudes are derived from the 
projection in Eqs. iJBJ and in Jacobi relative momenta. The ground state of the four-boson 
system has total angular momentum zero, therefore the spectator functions depends only on 
the modulus of the momenta. The integral equations are analogous to the ones shown in 
Ref. [26] (without three-body forces), but with the four-body Green's function substituted by 
their subtracted form as given in Eqs. © and Q. The numerical procedure followed Ref. [28]. 

The momentum scales in Eqs. © and (JJJ are a -1 , /i/3) and fxu) , consequently the tetramer 
binding energy for a -1 = depends on the momentum scales as 

S 4/M(3) = £i (M(4)/M(3)) • (8) 

In Table I, we present results for the tetramer ground state binding energy for \a\ — 00 and 
considering a fixed /i(3), such that we take fir^ as our momentum unit. From the numerical 
solution of Eqs. I@ and Q we obtain £4///^ for different values of /i(4)/M(3)- The trimer 
binding energy is fixed and given by the solution of Eq. ©, £3 = —0.00927^^ [23]. We 
clearly observe in this table the dependence of E4 on [lu) (for fixed /U(3)). As f-iu) / fJ,($) is 
varied from 1 to 20, the ratio E^jE^ changes from 5 to about 78 (considering two digits). 
In order to have £4 independent on the regularization scale fJ-u), one should require that for 
//(4)///(3) >> 1 the four-body scale would vanish as a physical scale. This is not observing in 
our results. 

In view of the scaling relation Eq. ©, the results presented in the table are valid for any 
/Lt (3 ). Using that E 3 = -0.00927/As in Eq. © we obtain a family of Tjon lines, where each 
one corresponds to a given value of H(4)/ From this result arises quite naturally the main 
claim of our work: the existence of a four-body scale, which is manifested through a family of 
Tjon-lines. The line slope is given by the ratio M(4)//*(3)- 

For equal three and four-body scales, i.e., nu) = /i(3) = 1, the tetramer binding energy is 
E 4 ~ 5 E 3 . The Tjon line [20] is E a = A.72(E t + 2.48) MeV (E a is the 4 He binding energy 
and Et the triton one), with the slope approached by our calculation. Also, it was recently 
found in Ref. [26] that E4 scales as ~ 5 E 3 , for equal bosons interacting through a zero-range 
two-body force, with a repulsive three-body potential needed to stabilize the trimer ground 
state energy against collapse. Our result for /j,^ — /i( 3 ) agrees with both Refs. [20] and [26]. 
The physical reason for this agreement is explained in the next two paragraphs. 

In the nuclear case, the nucleon-nucleon interaction is strongly repulsive at short distances 
and therefore the probability to have four nucleons simultaneously inside a volume ~ rj] is 
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Table I - Four-boson binding energy for E2 — and fixed /im ■ The three-boson energy is E3 = 
-0.00927/4,) [23], 



M(4)/M(3) 


1 


2 


4 


7 


10 


15 


20 




5.0 


7.8 


13 


22 


29 


51 


78 


E A /nf 3) 


-0.046 


-0.072 


-0.124 


-0.20 


-0.27 


-0.47 


-0.72 



quite small. For other potentials, like separable ones, parameterizing the nuclear force range, a 
similar result is found [29] , even without a short range repulsion due to the non- locality of the 
interaction. In these cases, presumably the four-nucleon scale itself has much less opportunity 
to be evidenced, i.e., the three- and four-body scales correspond to the same short-range 
physics, which in our calculation come from fiu\ = /*(3). 

In Ref . [26] , a three-body force is used to stabilize the shallowest three-body state, against 
the variation of the cut-off. Their results show that the Tjon line is reproduced when the three- 
body force is repulsive. Within this picture, the contribution of the three-body interaction in a 
configuration where the four bosons interact simultaneously is about four times more repulsive 
than having three-bosons in the force range and the other outside this range. Therefore, the 
dynamics of the short-range part of the interaction is manifested in the four-boson system 
mainly through small three-boson configurations, this should also qualitatively corresponds 
to the particular case jj,u) — /Lt/3) . 

Also, we should note that, in Ref. [30], considering the three-boson problem near a narrow 
Feshbach resonance, Petrov concludes that three-body observables depend only on the reso- 
nance width (which implies in an effective range of the order of A*^) and the scattering length. 
One may suspect that his conclusion can also have consequences to four and more particle 
observables. However, we would like to stress that, within our present coupled-channel for- 
malism, three and more body forces may appear, giving more freedom to the three and four 
atom systems. If one sticks to the Petrov's assumption, than the three and four body scales 
are equal (^(3) = /*(4)) and no independent behavior can be verified in practice. 

The general scaling form of s-wave three-boson observables, Eq. J2J) , can be generalized to 
four-bosons with the dynamics parameterized by the scattering length, trimer and tetramer 
binding energies. A s-wave four-boson observable will be strongly correlated to a, E3 and £4: 

04E,E 4 ,E 3 ,a) = \E^F 4 (J^ > ^ > ay/\E^j , (9) 

where O4 represents cither an observable at energy E or an excited tetramer energy (the 
dependence on E does not appear in the last case) . The exponent 77 gives the correct dimension 
to O4. There is the possibility of a new renormalization group limit cycle due to the new scale, 
and for /im — * 00 one may speculate that the correlation given by Eq. would approach a 
universal function. 

In summary, by studying a four-boson system with a zero-range interaction, we conclude 
that a four-body scale is necessary when the two and three-boson scales are fixed. Extrapolat- 
ing this result to more particles, we suggest that one new experimental information is required 
for each new particle added to the system near a Feshbach resonance, to include effects from 
many-body forces from the one-channel reduction. In the particular case, when only two-body 
forces are present, all scales are equal and can be for example associated with the effective 
range. However, as more accuracy is needed to describe the physics near a Feshbach resonance 
a distinction between 3 and 4 body scales must be introduced. 
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